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Abstract 

The scalar susceptibility (x s ) of QCD, which represents the response of the chiral condensate to a 
small perturbation of explicit chiral-symmetry breaking (m ^ 0), is investigated within the nonlocal 
chiral quark model (NL%QM) based on the instanton vacuum configuration for Nf = 2. We also 
take into account 1/N C meson-loop (ML) corrections including scalar and pseudoscalar mesons. It 
turns out that the chiral condensate is modified to a large extend by the ML corrections in the 
vicinity of m = 0, whereas its effect becomes weak beyond m ~ 100 MeV. As numerical results, we 
t— 5 find that Xs = — 0.34GeV 2 with the ML corrections and 0.18 GeV 2 without it, respectively. From 

these observations, we conclude that the ML corrections play an important role in the presence of 
finite current-quark mass. 
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I. INTRODUCTION 



Among various QCD susceptibilities, which associate with the response of the QCD to 
external sources, scalar (chiral) susceptibility, Xs plays an important role in understanding 
the pattern of explicit chiral-symmetry breaking. One can define it as follows: 

x. = E m) 



Sin 

flavor 



(i) 

m=0 



where (qq) and m stand for the chiral condensate and current-quark mass, respectively. 
Note that this nonperturbative QCD quantity is also of great importance in various QCD 
environments, such as the vacuum, finite p and/or T, and quark-gluon plasma, since it 
is deeply related to the phase structure of QCD with the physical quark mass. To date, 
there has been a good amount of related works, for example, QCD (LQCD) simulation pQ, 
Schwinger-Dyson method [21 [3], multiflavor Schwinger model [H [5], linear-a model [61 [7], 
Nambu-Jona-Lasinio (NJL) [HI9], perturbative QCD (pQCD) [ID] , and so on. 

In the present work, we would like to investigate the scalar susceptibility Xs in the vacuum 
by employing an effective chiral action, derived from the instanton vacuum configuration, 
i.e., the nonlocal chiral quark model (NLxQM). This highly constrained framework is char- 
acterized by inter-instanton distance (R ~ lfm) and average size of instantons (p ~ 1/3 fm), 
which provides a scale of the model, A w 1/p. In many practical applications of the model, 
large N c and chiral limit have been taken into account, which simplify the usage of the model 
considerably and has given good descriptions of light-flavor sector [Tl| fl~2| H3]. Attempts to 
overcome these two limiting cases (N c — > oo and m — > 0) have been done already extensively 
in Refs. [HI [151 EH EH (HI EE 120], in which the finite current-quark mass was explicitly 
treated [TU [T31 HUl [EH HE] , and meson-loop (ML) corrections, showed a sizable modification 
to the leading N c contribution [TDJ ED]. From these previous works, it also turned out that 
the ML corrections are crucial in treating the finite current-quark mass. Hence, it is natural 
and critical for us to consider the ML corrections carefully in the present work, since we are 
interested in the response to a finite value of m. 

We first compute the m-dependent constituent-quark mass, M (m), by solving saddle- 
point equations with the ML corrections, scalar and pseudoscalar mesons, generated by 
bosonizing an instanton- induced effective action for Nf = 2. We observe that the ML 
corrections gives considerable modification in M for the region below m ~ 100 MeV, in 
comparison to that from the leading-order contribution only. In addition, for computing Xs, 
chiral condensate is also computed, resulting in {iq\ 'q) m =i50MeV / '{^ Q) m=o ~ 0.7 with the 
ML corrections, whereas the ratio becomes about 0.4 without it. Note that the empirical 
value for (is^ s) / (ivJu) ~ (is* s) / (id) d) 0.8 [2T| . According to Eq (JT1), in Minkowski space, 
we reaches Xs — — 0.34GeV 2 with the ML corrections and 0.18GeV 2 without it. From 
these observations, we can conclude that the ML corrections plays an important role in 
considering small but a non-zero current-quark mass, showing obvious change of the scalar 
susceptibility in its sign and magnitude. As for the strength of Xs, since scalar mesons may 
not be sensitive to a small perturbation of the current-quark mass, the computed values for 
Xs are reasonable in comparison to pseudoscalar susceptibility ~ 1 GeV 2 [B] . 

The present work is organized as follows: In Section II, we briefly review the general 
formalism of the NL^QM with the ML corrections. In Section III, we present numerical 
results and discussions for the chiral condensate and scalar susceptibility. The final Section 
is devoted to summary and conclusion. 
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II. EFFECTIVE CHIRAL ACTION FROM THE INSTANTON VACUUM 



In this Section, we would like to provide a brief introduction of the present theoretical 
framework. As mentioned in the previous section, we utilize the nonlocal chiral quark model 
(NLxQM), based on the instanton vacuum configuration. Note that all calculations will be 
performed in Euclidean space hereafter. Otherwise, we will mention. First, we start with a 
Dirac equation for a quark in the presence of instanton background fields: 

W + im-fldVf^WtW, (2) 

where m designates the current-quark mass. Ajj stands for a singular- (anti) instanton solu- 
tion as follows: 



2rf??x v 

AlI » {x) ~ x*(xi + p*y 



Here, rj® u indicates the 't Hooft symbol as given in Ref. [22], whereas p the average size 
of (anti)instantons. I and I stand for instanton and anti-instanton contributions, respec- 
tively. If we assume that low-energy nonperturbative QCD properties are dominated by the 
following fermion (quark) zero-mode (Ao = generically) , 

[ip + im- Ail] */? = 0, (4) 
we can write a quark propagator with the zero- mode solution, ty^j , approximately as 

1 tfWtfg? 
ip + im — /Ijj vm 

where S denotes a free-quark propagator: 

So = — ( 6 ) 
ip + im 

Eq. (|5]) can be written alternatively via the Fourier transform of the zero-mode solution: 

g — (7) 

ip + im + iM{id,m) } 

where momentum and current-quark mass dependent quark mass, M(id, m) is defined by 

M(id, m) = M (m)F 2 (id). (8) 

Here, Mq means constituent-quark mass. Moreover, since we are interested in finite current- 
quark mass, Mo has been taken into account as a function of m implicitly. The current-quark 
mass dependence of M will be discussed in detail later. The form factor F, which comes 
from the Fourier transform of the quark zero-mode solution, reads: 



F(t) = It 



mK.it) -h{t)K {t)- hit)Kl{t) 
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where K n and I n are modified Bessel functions. With all ingredients having considered, one 
can construct an effective partition function, which brings about the quark propagator given 
in Eq. Q, as follows [22]: 



Zefi = dX Dip Dip^ exp 



d A x^(ip + im)i) + X(Y + + Y~) + N (in 



N 



XV M 



- 1 



(10) 



where the A designates the Lagrangian multiplier to make the action exponent [TH]. More- 
over, we have introduced an arbitrary massive parameter M. to make the argument of the 
logarithm dimensionless and considered the numbers of the instanton and anti-instanton are 
the same, Nj = Nj = N. Since we are interested in Nf = 2, the 't Hooft interaction can be 
written as: 



Y 2 ± = J d 4 z dU J d 4 xi d A X2 d A x 3 d A x^ 



x 

X 



^{x-y) L)R {id + imf^Uxi - z)^\x 2 - z)(id + im) ± ip Lfi {x 2 ) 



n 

^(x 3 ) LtR (id + im) T ^(x 3 - z)^\x 4 - z)(id + im) ± ^L,K(^4) 



n 

(o)- 
n 



(11) 



Here, the subscripts L and R represent the chirality of quarks and z the instanton center. 
We introduce a notation for simplicity as x ± = x M (±zcr, I2X2) [23] . 

An effective action can be obtained from the effective-partition function defined in 



Eq. (10). Note that the 't Hooft interaction given in Eq. (11) can be casted into four-quark 
instanton-induced interactions by integrating over color orientation U. From this Nambu- 
Jona-Lasinio (NJL) model-like four-quark interactions, one can construct an effective chiral 
action as a functional of quarks, and scalar and pseudoscalar mesons, $ = ($1, $2, $3, $4) = 
(a, 77, er, <j>), by an exact bosonization process for Nf = 2 [19, 22J. The meson fields are nor- 
malized as $ 2 = a 2 + f] 2 + cr 2 + 2 . Since there must be meson fluctuations around their 
classical paths, one can consider an additional contribution to the effective chiral action from 
these fluctuations as meson-loop (ML) corrections. Practically, we can write the effective 
chiral action with the ML corrections as follows: 

S eff = S e L ff ° + S e N ff LO , (12) 
where the leading-order (LO) term is obtained by bosonizing the effective-partition function 



of Eq. ( 10 ) [in] , resulting in a linear-a type action: 

-Nln-£— + 1 + 2 fd A x (a 2 + ? 1 2 + <T 2 + cj> : 
XV Jvl J ^ 



°eff 



— Spin 



+ im + iM F(id) (a + 7577 + it ■ cr + ry 5 r • 0) F(id) 
i0 + im 



(13) 



where Sp denotes a functional trace, / d 4 x Tr c j n (x\ ■ ■ ■ \x), in which the subscripts c, /, and 
7 indicate color, flavor, and Dirac-spin indices, respectively. Note that satisfies LO 
saddle-point equations, dS^/d^i = for i = 1 ~ 4. M is defined by \/X(2np) 2 /2N C in the 
large N c limit and equals to M /$i [P9"] . 

Now, we are in a position to consider the meson fluctuation explicitly, $ — > $ + $', which 
corresponds to l/N c corrections as pointed out in Refs [T9l EO] . This next-to-leading-order 
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(NLO) contribution can be derived straightforwardly by functional derivatives of S^s with 
respect to the meson fields: 



qNLO 

°eS,i 



-Sp $ In 



= -Spin 
2 F 



5<!> i (x)5<S> j (y)_ 
-Tr 



M F 2 (id) 



M F 2 (id) 



Ti 

_ip + im + iM F 2 (id) T ■ <3> ip + im + iM F 2 (id) T • $ 



-r, : 



(14) 



where Sp$ runs over the meson fields. Also we define T = (1, 75, ir, i^r) and Tr = Tr c j i7 . 
Hereafter, we will take into account the finite fluctuation only in the direction of the isoscalar 
a field, corresponding to the spontaneous breakdown of chiral symmetry (SB%S, (a) 7^ 0), 
whereas other fluctuations are set to be zero, resulting in T ■ $ — > a. By combining the LO 
and NLO contributions, we finally arrive at the following expression for the effective chiral 
action with the meson-loop (ML) corrections: 



oLO+NLO _ Ari N 

S ° s ~ ~ Nln XVM 



+ 1 + 2/ dW-Sp cJi7 ln 



ip + im + iM F 2 (id) 



+ im 



LO 



i=i z 



4- ^Tr 



M F 2 (id) M{id) 



ip + iM(id) l ip + iM(id) 



(15) 



NLO 



Here, M(k) stands for m + M (m)F 2 (k) = m + M(k). It is more convenient to define the 
effective chiral action in momentum space in order to compute physical quantities: 



oLO+NLO 



C + AHnA + 2 / d 4 xa 2 -V [ y^Trln 

J J (27T) 4 



i=l Z 



(27^ 



In 



1 - 



(27T) 

d 4 A; 



4(72 7 (2tt) 



Tr 



A4 r 

Pa Pb 



(16) 



where p a ,6 = $a,& + iM a j and i a = k* + im with k a = k and kb = k + q. The first term, 
C represents irrelevant constant terms. By differentiating the effective chiral action with 
respect to A and <x, one is led to two-individual saddle-point equations, which satisfy the 
following vacuum equations 



0£LO+NLO 



^jjLO+NLO 

da 



= 0, 



(17) 



N 
V 



LO 
/1 7 



4 7 (2tt) 4 W 4^y (2tt) 4 



LO 



NLO 

d 4 fc 



^ 2 -I^G t (k,q) 



(27T) 



(18) 



(19) 



NLO 
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In the above equations, we have employed the following notations for simplicity: 



F(k) = Tr 



iM a 

Pa 



G t (k,q) = -Tr 



Pa 'P, 1 



Hi(k,q) = 7 Ti 



MaMa Mb 
Pa Pa . *A % 



(20) 



Detailed evaluations of F, Gi, and Hi are given in Appendix. By incorporating Eqs. (18) 
and (19), M (m) can be calculated numerically. However, instead of doing that, considering 
the fact that the LO term in Eq. (19) is lager than the NLO term by a factor 10 ~ 20, we 
approximate a 2 in the denominator of Eqs. ( 18 ) and ( 19 ) as 

1 



a 



a 2 



d A k 



4 J (2tt) 



F(k). 



(21) 



We note that this replacement also has to do with a proper LO chiral behavior of a meson 
propagator at small meson momentum (q — > 0), that is, the Gell-Mann-Oakes-Renner (GOR) 
relation, as indicated in Refs. [TTJ [21] : 



11;^) = 4 



a 



d 4 k 



Gi(k,q) 



d A k 



[F(k)-AGi(k,q)] 



(22) 



(2vr) 4 

As for the pion (i = 4) for Nf = 2 and « 0.14 GeV, from Eq. (22), we obtain the 
following relation in the leading order 0(m) in Euclidean space: 



n; x (o) 



. m N c N f 



d A k 



M a 



k 2 + Ml 



GOR relation 

where we take f w ~ 93 MeV, and the pion propagator reads 
„ M 1 1 



V 2 [q 2 + ml 



(23) 



(24) 



Although we have used an analytical LO expression for (iq^q) in Eq. (23) without expla- 
nation, it will be discussed in detail soon. Plugging a 2 into Eq. (18) and using instanton- 

(0.2 GeV) 4 



packing fraction N/V ~ (0.2 GeV) [TT], one can have M as a function of m. In practical 
calculations, we make use of the following parameterized form factors, which mimics that in 
Eq. ^ qualitatively well, because of numerical simplicity: 



M(k) = M 



2A 2 



|fc|<2.0GeV, M(k) = M 



6 



(\kW 



\k\ > 2.0 GeV, 



(25) 



2A 2 + k 2 

where we choose A ps 1/p ps 0.6 GeV. 

In the left panel of Fig. [TJ we draw numerical results for Mq as a function of m for the 
LO (without ML corrections) and LO+NLO (with ML correction) contributions, separately. 
It turns out that M LO = 253 MeV and M LO+NLO = 179 MeV at m = 0. Being distinctive 
from M LO , which decreases monotonically with respect to m, M LO+NLO increases and then 
decreases, indicating the different portion of vacuum effects depending on m. In the right 
panel of Fig. [TJ we depict a 2 in the same manner with M . As for the LO contribution, 
we have (a LO ) 2 ps N/(2V) « 8.0 x 10~ 4 GeV 4 . On the contrary, (^lo+nlo^ starts from 
about 3.0 x 10" 4 GeV 4 and come closer to (cr LO ) 2 as m grows. We note that the present 
results for Mq and a 2 with ML corrections are in qualitative agreement with those given in 
Refs. [T9l 1201 124"] . employing the same theoretical framework. In their works, however, Mq 
and a 2 were computed starting from the chiral limit (m = 0), and the chiral corrections was 
evaluated. On the contrary, we obtain them directly from the saddle-point equations with 
the fairly good approximation, a 2 — > b 2 in the meson propagator, 11$.. 
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400 



9.0 




100 50 100 150 200 1 '°0 50 100 150 200 

m [MeV] m [MeV] 

FIG. 1: Mq (left) and a 2 (right) as functions of m. We draw it for the LO (dashed line) and 
LO+NLO (solid line) contributions, separately. 



III. SCALAR SUSCEPTIBILITY WITH MESON-LOOP CORRECTIONS 



In this section, we will calculate chiral the condensate and scalar susceptibility utilizing 
the effective chiral action with the ML corrections, derived in Section II. Here, we write the 
relevant terms of the effective chiral action: 



qLO+NLO 

° e fr 



-V 



( -, — rrTrln 


\Pa\ 


/ (2tt) 4 


. 4a _ 



+ E 



V f d 4 q 



i=i 



(27T) 



hi 



1 - 



1 

4^ 



d 4 k 



Tr 



U, Mb 

Pa Pb 



(26) 



'cS 



where d a = ]t a + im. The chiral condensate can be derived straightforwardly from ^lo-nlo 
by a functional derivative with respect to m, resulting in: 

1 a<7 L P+NLO 

V iV f UIIL 

J^H l (k,q)/M a 



VN f 
-AN r 



dm 
d A k 



(2vr) 



Ma 



m 



+ 



N 



E 

/ i=l 



d^q 



* 2 -fjSr*Gi(k,q) 



(27) 



LO 



NLO 



where D\ = D a D\ and d\ = d a d\. Note that a 2 in the denominator of the NLO term 
has been replaced by a 2 as done previously. In Fig. [2J we draw the numerical results of 
(iq^q) for the LO and LO+NLO contributions, separately. At m = 0, it turns out that 
(zgtg)LO = (253 MeV) 3 , whereas (iq^q) LO+NLO = (210 MeV) 3 , showing about 20% difference. 
Interestingly, (ig"l"g) LO+NLO is saturated once around m ~ 20 MeV, then decreases, being 
different from the LO contribution, which decreases monotonically. Two contributions start 
to behave in a similar manner beyond m ~ 100 MeV as Mq does. 

Although we have been working with Nf = 2, it is meaningful to examine the chiral 
condensate away from m — 0, to see the effects of the ML corrections for explicit chiral 
symmetry breaking. Empirically, the ratio between the condensates of light and strange 
quarks was estimated to be (ss) / (uu) ~ (ss) / (dd) = 0.8±0.3 [25] and 0.75±0.12 [26]. This 
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FIG. 2: (iq^q) as a function of m. We draw it for the LO (dashed line) and LO+NLO (solid line) 
contributions, separately. 



ratio in the present framework at m = and m = 150 MeV, being computed, it turns out to 
be 0.71, which shows relatively good agreement with the empirical values. In contrast, the 
LO contribution gives the ratio ~ 0.43. From these observations, we can conclude that the 
ML corrections improve the model noticeably in the presence of finite current-quark mass. 

Now, we are in a position to compute the scalar susceptibility Xs, incorporating Eqs. ([I]) 
and (27). In Minkowski space, Xs reads: 



Xs 



E 

flavor 



E 



dm 



d 4 q 



-N; 



d{iq^q) 



771=0 



dm 



AN r N 



f 



d 4 k 



1 dHi(k,q) 
M a dm 



7Tt=0 



777 =0 



d 4 k 



D 2 



2M? 



n4 



1 



d 4 k 
(2tt) 4 



LO 

Hi(h,q) dGj(k,q) 
M a dm 



777 =0 



a? 



/ (2tt) 4 ^'(^) l) 



(28) 



where p a0 = fc a + iM a and erg 
m are evaluated as: 



NLO 



a 2 at m 



0. The derivatives of F and G with respect to 



dHi(k, q) 



dm 
dGi(k,q) 



Tr 



dm 



777 =0 



771=0 



M a M a 



M h 



PaO PaO PaO P 



hi) 



-Tr 



Mh 



r,;^r, ; 



PaO PaO P, 



hi) 



(29) 



As for the numerical results, we obtain x^° ~ 0.18 GeV 2 and Xs° +NL ° ~ -0.34 GeV 2 in 
Minkowski space, showing substantial modification from the ML corrections. Especially, 
the sign difference is crucial due to the obviously different behaviors of M near m = 
as shown in the left panel of Fig [TJ Note that Xs can be estimated from the results of 
Ref. [24J, giving about —0.50 GeV 2 , which is compatible with ours. Here, an explanation for 
the strength of Xs is given; since scalar mesons must be insensitive to a small perturbation 
of the current-quark mass being considered, their relatively heavy masses, the computed 
values for x are reasonable in comparison to the pseudoscalar susceptibility, estimated as 
~ 1 GeV 2 in Ref. [6] . All the numerical results up to now are listed in Table. [| 
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M [MeV] 


a 2 [GeV 4 ] 


\{iq*q)\s [MeV] 


{{iqUyfi [MeV] 




Xs [GeV 2 ] 


LO 


353 


8.00 x 10" 4 


255 


192 


0.43 


0.18 


LO+NLO 


179 


2.88 x 10" 4 


210 


187 


0.71 


-0.34 



TABLE I: Various physical quantities computed at m = 0. Here, we take m = 150 MeV for (iq^q)'. 

IV. SUMMARY AND CONCLUSION 

In the present work, we have investigated the scalar susceptibility of the QCD vacuum, 
X s and related nonperturbative quantities, such as the chiral condensate (iq^q). For this 
purpose, we employed an effective chiral action, derived from the instanton vacuum configu- 
ration in Euclidean space. Since it has been known that the meson-loop (ML) contribution, 
corresponding to the large N c corrections, play important roles in a system beyond the chi- 
ral limit and, since we were interested in the response of the chiral condensate to a small 
perturbation of current-quark mass, we added the ML corrections as a NLO contribution 
to the LO effective chiral action. By virtue of field-theoretical functional methods, the ML 
corrections were derived from the LO contribution by fixing meson fluctuations along the 
direction of the isoscalar o field. 

The saddle-point equations were obtained by differentiating the effective chiral action 
with respect to the external parameters, A and a, giving current-quark mass dependent 
constituent quark mass M , self-consistently. Owing to the ML corrections, M behaves very 
differently in the vicinity of m = from that without the corrections. This difference effected 
also on the chiral condensate as well as the scalar susceptibility, substantially. Especially, 
we observed in Minkowski space that Xs = 0.18 GeV 2 without the ML corrections, whereas 
Xs = —0.34 GeV 2 with it, showing the critical sign difference. We note that this small 
value, in comparison to that for pseudoscalar susceptibility ~ 1 GeV 2 , is moderate, taking 
into account that heavier scalar mesons might be insensitive to a small perturbation of 
current-quark mass. 

Consequently, we can verify that the ML corrections are critical in considering a system 
beyond chiral limit as shown in the drastic change in scalar susceptibility for instance. As a 
perspective, we are also concerned in the extension of the present results to a system with 
Nf = 3, and finite density and/or temperature. These extensions must be important in 
exploring QCD phase structure and finding QCD critical point as exemplary topics. 
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Appendix 

The functions F(k), Gi(k, q), and H^k, q) are evaluated for i = 1 ~ 4: 

F(h , AN c NfM a M a 
{) Dl ' 

N c N f M a M b (k a -k b -M a M b ) 



DIDl 

-N c N f M a M b {k a ■ k b + M a M b ) 
DIDl 

-3N c N f M a M b (k a ■ k b - M a M b ) 
n2n2 ' 

U a U h 

3N c N f M a M b (k a ■ k b + M a M b ) 
DIDl 

N c N f M 2 M b [M b k 2 a + 2M a (k a ■ k b ) - M 2 a M b ] 

U a U b 

N c N f MlM b [M b kl - 2M a (k a ■ k b ) - M 2 a M b ] 

-3N c NfM 2 a M b [M b k 2 a + 2M a {k a ■ k b ) - M 2 a M b 
DiD 2 

-3N c N f M 2 M b [M b k 2 a - 2M a {k a ■ k b ) - M 2 M b 
DiD 2 
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